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SAMENVATTING 

Scheepsbewegingen  kunnen  zeeziekte  veroorzaken.  Het  zijn  daarbij  vooral  de  verticale 
bewegingen  met  frequenties  onder  de  1  Hz  die  een  provocatief  karakter  hebben.  In  dit 
rapport  wordt  een  samenvatting  en  de  beschrijving  van  een  computer-implementatie  gegeven 
van  een  model  zoals  geformuleerd  door  McCauley  e.a.  (1976)  en  een  model  volgens  Griffin 
(1990),  waarmee  het  percentage  zeezieken  berekend  kan  worden  op  basis  van  een  gegeven 
scheepsbewegingsprofiel . 
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Ship  motions  and  “Motion  Sickness  Incidence”  (MSI):  Estimations  of  the  percentage 
sick  people  based  on  vertical  motion 

J.E.  Bos 


SUMMARY 

Ship  motions  can  induce  seasickness.  Especially  the  vertical  movements  with  frequencies 
below  1  Hz  are  very  provocative.  This  report  gives  a  summary  and  a  description  of  a 
computer-implementation  of  a  model  formulated  by  McCauley  et  al.  (1976),  and  a  model  by 
Griffin  (1990),  to  predict  the  percentage  of  sick  people  given  a  motion  time  history. 
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1  INLEIDING 

Bij  het  normale  gaan  en  staan  kan  de  (gezonde)  mens  zich  probleemloos  ten  opzichte  van  de 
zwaartekracht  orienteren  en  zijn  evenwicht  bewaren.  Begeeft  hij  zich  op  ander  terrein, 
bijvoorbeeld  op  zee,  en  ondergaat  hij  andere  bewegingen  dan  normaal  op  de  wal,  dan 
kunnen  problemen  ontstaan.  Deze  problemen  varieren  van  een  gestoord  houdingsevenwicht 
tot  onwelbevinden,  eventueel  gepaard  gaande  met  braken  (zeeziekte).  De  mate  van  het 
onwelbevinden  kan  sterk  per  persoon  verschillen.  Het  percentage  mensen  dat  moet  braken 
wordt  de  “Motion  Sickness  Incidence”  (MSI)'  genoemd.  Het  is  bekend  dat  bij  verticale 
bewegingen  mensen  meer  ziek  worden  dan  bij  andere  bewegingen.  Om  hiermee  bij  scheeps- 
ontwerp  rekening  te  kunnen  houden  is  in  het  verleden  getracht  dit  percentage  zeezieken  te 
kennen  als  functie  van  een  gegeven  bewegingsprofiel  (zie  bijvoorbeeld  Wertheim  e.a., 
1995a).  Naarmate  de  trend  doorzet  dat  schepen  met  steeds  minder  bemanning  worden 
uitgerust,  wordt  het  steeds  belangrijker  dat  het  functioneren  van  bemanningsleden  zo  weinig 
mogelijk  gehinderd  wordt  door  de  scheepsbewegingen.  De  beschrijving  van  de  MSI  als 
functie  van  de  beweging  is  het  onderwerp  van  dit  rapport.  De  twee  in  de  literatuur  meest 
genoemde  modellen  zullen  hier  behandeld  worden.  Tevens  is  van  deze  modellen  bij 
TNO-TM  een  computer-programma  gemaakt. 

Het  eerste  model  komt  van  McCauley  e.a.  (1976).  Zij  hebben  in  een  scheepsbewegingssimu- 
lator  619  (mannelijke)  proelpersonen  blootgesteld  aan  verticale  sinusoidale  bewegingen  met 
RMS-amplitudes  varierend  van  0.0278  tot  0.555  G  en  frequenties  van  0.083  tot  0.7  Hz,  met 
een  duur  tot  2  uur.  Gekeken  werd  welk  percentage  van  deze  proelpersonen  binnen  die  tijd 
moest  braken.  Deze  M57-scores  bepalen  de  parameters  van  een  volgens  hen  plausibel  model. 
Dit  model  heeft  slechts  drie  input-parameters:  de  amplitude,  de  frequentie  en  de  duur  van  de 
beweging.  Het  McCauley-model  ligt  ten  grondslag  aan  de  ISO  2631/3  norm. 

Het  tweede  model  is  gebaseerd  op  het  werk  van  Griffin  (1990)  volgens  de  “British 
Standard”  6841.  Dit  werk  beschrijft  het  welbevinden  van  de  mens  onder  allerlei  vormen  van 
vibratie.  Zeeziekte  kan  dan  opgevat  worden  als  een  extreme  vorm  van  onwelbevinden  bij 
vibraties  met  een  lage  frequentie  (0.1— 0.5  Hz).  Dit  model  houdt,  in  tegenstelling  tot  het 
model  van  McCauley  e.a.,  expliciet  rekening  met  de  mogelijkheid  dat  een  beweging  uit 
meerdere  frequenties  bestaat.  Net  zoals  bij  McCauley,  op  wiens  werk  een  deel  van  zijn 
model  gebaseerd  is,  verschilt  de  zeeziektegevoeligheid  per  frequentie.  Griffin  gebruikt 
daarnaast  data  over  het  welbevinden  van  passagiers  op  veerboten.  Verder  wordt  zijn  model 
verondersteld  geldig  te  zijn  voor  bewegingen  met  een  duur  tot  6  uur.  De  input  van  het 
Griffin-model  is  dus  een  compleet  verticaal  bewegingsprofiel,  d.w.z.  voor  iedere  t  in  het  te 
onderzoeken  tijdsinterval  is  de  amplitude  van  de  verticale  beweging  gegeven. 

Voor  de  volledigheid  zij  nog  verwezen  naar  een  rapport  van  de  Technische  Hogeschool 
Delft  (Vermeij,  1985),  over  onderzoek  en  voorspelling  van  zeeziekte.  Vermeij  komt  in  dat 
rapport  tot  de  conclusie  dat  vooral  een  model  van  Krappinger  (1983)  en  een  van  Smith 
(1982)  ter  voorspelling  van  de  MSI  mogelijk  bruikbaar  zijn.  Het  model  van  Krappinger  is 
echter  op  de  parameterkeuze  na  identiek  aan  het  model  van  McCauley.  Het  model  van  Smith 


Vanwege  hun  algemeen  gebruik  wordt  in  dit  rapport  een  aantal  Engelse  naamgevingen  gehandhaafd. 
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overlapt  dat  van  Griffin  waaraan  toegevoegd  de  gevolgen  van  mogelijke  interacties  tussen 
verschillende  frequentie-componenten.  Omdat  er  niet  voldoende  experimentele  gegevens 
voorhanden  zijn  om  de  parameters  te  schatten  die  deze  interacties  kwantificeren,  is  deze 
toevoeging  op  dit  moment  echter  niet  bruikbaar. 


2  MODEL  VOLGENS  McCAULEY 
2.1  Theorie 

Algemeen  wordt  aanvaard  dat  het  voornamelijk  de  verticale  bewegingen  met  frequenties 
tussen  ongeveer  0.1  en  0.7  Hz  zijn,  die  zeeziekte  veroorzaken.  Verticaal  wil  bier  zeggen 
parallel  aan  de  zwaartekrachtrichting.  Dit  wil  dus  zeggen  dat  alle  andere  scheepsbewegingen 
(de  overige  twee  translaties,  de  drie  rotaties  en  dompbewegingen  trager  dan  0.1  Hz  en 
sneller  dan  0.7  Hz)  verondersteld  worden  in  het  geheel  geen  zeeziekte  te  veroorzaken. 
Zowel  McCauley  als  Griffin  modelleren  dan  ook  de  MSI  enkel  als  functie  van  deze  verticale 
bewegingen. 


time  (min) 

Fig.  1  MSI  als  fimctie  van  de  tijd  bij  drie  verschillende  G-waarden  (RMS),  en 
66n  vaste  frequentie  (f=0.25  Hz).  Getrokken  lijnen  zijn  modelberekeningen. 


Op  basis  van  alle  gegevens  (het  aantal  spugers  per  amplitude  per  frequentie)  verzameld  door 
McCauley  e.a.  in  de  bewegingssimulator  van  het  Office  of  Naval  Research/Human  Factors 
Research,  Incorporated,  Goleta,  California,  kwamen  zij  tot  de  volgende  conclusies.  Het 
bleek  dat  de  zeeziektegevoeligheid  als  ftinctie  van  de  logaritme  van  de  tijd  eerst  langzaam 
toeneemt,  met  een  maximale  toename  bij  een  zekere  (log  van  de)  tijd  n„  waarna  voor 
grotere  t  altijd  een  verzadiging  optreedt  (zie  Fig.  1). 
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Fig.  2  MSI  (%)  versus  frequentie  en  stimulus  amplitude.  Bolletjes  geven  de 
meetpunten  weer,  lijnen  de  optimale  fit  door  deze  punten  (zie  tekst). 


Een  soortgelijke  afhankelijkheid  werd  waargenomen  voor  de  aangeboden  versnell ingen,  zij 
het  dat  hierbij  de  versnelling  waarbij  de  maximale  toename  optreedt  afhankelijk  bleek  te  zijn 
van  de  stimulusfrequentie  (/).  Deze  bevindingen  staan  grafisch  weergegeven  in  Fig.  2.  De 
MSI  is  bier  steeds  het  percentage  van  een  groep  dat  binnen  de  gegeven  duur  van  twee  uur 
meet  overgeven,  dan  wel  de  kans  dat  een  individu  dit  doet.  Nadrukkelijk  zij  vermeld  dat  de 
MSI  hier  dus  een  cumulatieve  maat  is,  wat  wil  zeggen  dat  de  MSI  het  totaal  aantal  spugers  in 
die  gegeven  tijd  geeft.  Om  de  in  Fig.  1  en  2  weergegeven  data  in  66n  mathematisch  model 
te  vatten  worden  eerst  de  tijds-  en  versnellingsafhankelijkheid  benaderd  met  een  cumulatieve 
normale  verdeling^: 

MSI  ~  ^(z)  =  -Lfe  ^  dz  (1) 

27r  j 


met  z  =  z,  voor  de  tijdsafhankelijkheid  en  z  =  voor  de  versnellingsafhankelijkheid: 


_  log/  -  p, 


en 


logfl  -  p^(/) 
P^a 


(2) 


^  Citaat  uit  McCauley  e.a.  (1976,  p.  43):  “The  genesis  of  his  model  has  no  doubt  been  influenced  by 
the  fondness  of  one  of  the  authors  for  the  concept  of  logarithmic  responses  to  physical  stimuli,  and  by 
his  trust  in  the  Central  Limit  Theorem  to  cause  normal  distributions  to  appear  conveniently  when 
needed”. 
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Hierin  stelt  p  een  constante  voor  zodanig  dat  de  waargenomen  MSI  gelijk  is  aan  de  voor- 
spelde.  Deze  p  wordt  door  McCauley  e.a.  gegeven  door  p  =  V'(l-p^)  waarin  p  de  correlatie 
tussen  waargenomen  en  voorspelde  (log  van  de)  versnelling  is  waarbij  gebraakt  wordt.  De 
logaritme  is  verder  die  met  bet  grondtal  10.  De  sigma’s  geven  de  breedte  aan  van  bet  gebied 
waarin  de  Af5/-toename  bet  grootst  is.  Samengevat  leidt  dit  tot: 


MSI  =  1004>(a,0  =  - — - 


De  afhankelijkheid  van  de  versnelling  waarbij  de  MSI  maximaal  toeneemt  als  functie  van  de 
frequentie  is  door  McCauley  e.a.  als  volgt  nader  gespecificeerd: 

De  in  (3)  gegeven  bivariate  verdeling  is  een  lastige  functie  (niet  te  vinden  in  standaard 
tabeilen  bijvoorbeeld).  Wanneer  we  echter  aannemen  dat  en  z,  statistisch  onafhankelijk  van 
elkaar  zijn,  dan  mag  de  dubbele  integraai  in  (3)  gesplitst  worden  naar  bet  produkt  van  de 
twee  afzonderlijke  integralen,  waarbij  de  substitutie 


vT^ 


bet  probleem  met  de  kruisterm  in  de  exponent  van  vergelijking  3  opiost.  Samengevat  geeft 
dit  nu: 


M5/ =  100  4)(z„)$(z/)  (6) 

Met  bebulp  van  de  bulk  aan  data  en  bet  bovengegeven  model  zijn  door  McCauley  e.a.  de 
parameters  k„  k2  en  k;  voor  de  bepaling  van  a^,  ij.„  a,  en  p  geschat.  Het  resultaat  staat  in 
Tabel  I  gegeven.  Met  deze  parameters  wijkt  bet  model  niet  meer  dan  3.8%  af  van  de 
werkelijk  gemeten  data  (sd=:2.4%). 


Tabel  I  McCauley  model-parameters. 


p.  =  0.87  +4.36  log/+2.73  logV 

fii  =  1.46 

a,  =  0.47 

(T,  =  0.76 

p  =  -0.75 

2.2  Model  implementatie 

Met  bebulp  van  vergelijkingen  2  en  4-7  en  de  bulp  van  een  tabellenboek  om  4>(z) 
[=  l-$(-z)]  te  bepalen,  kan  de  MSI  (bandmatig)  berekend  worden.  Voor  de  functie  ^ 
bestaan  ecbter  goede  benaderingen  die  met  een  computer  eenvoudig  zijn  te  berekenen.  Dus 
kan  een  eenvoudig  computerprogramma  een  gewenst  alternatief  zijn  voor  de  bandmatige 
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berekening  met  gebruikmaking  van  een  tabellenboek.  McCauley  e.a.  geven  in  bun  rapport 
een  voorbeeld  van  zo’n  programma,  welke  met  de  toevoeging  van  de  in-  en  output-voorzie- 
ning  in  FORTRAN  voor  een  MS-DOS  machine  door  TNO-TM  gerealiseerd  is.  Een 
voorbeeld  van  de  vragen  en  resultaten  zoals  dit  programma  die  stelt  en  weergeeft  staat  in 
Bijlage  A.  Met  nadruk  zij  gewezen  op  dat  waar  Griffin  (1990)  Sl-eenheden  gebruikt, 
McCauley  e.a.  dat  niet  doen.  In  de  geimplementeerde  versie  van  het  McCauley-model  pas  ik 
me  hierbij  aan,  en  dient  de  versnelling  in  G-eenheden  te  worden  opgegeven,  de  tijd  in 
minuten  en  de  frequentie  gewoon  in  Hz. 

2.3  Resultaten 


Met  een  aangepaste  versie  van  bovenbeschreven  programma  is  voor  ieder  interval  van  0.1 
<  F  <  0.5  Hz,  0.3  <  A  ^  5  m/s^  en  20  ^  T  <  120  minuten  de  MSI  berekend  bij  drie 
waarden  voor  de  overige  twee  variabelen.  Deze  resultaten  staan  weergegeven  in  Fig.  3. 


^  (Mz)  A  (•/a2)  T  (nin) 

Fig.  3  MSI  versus  frequentie  (F),  amplitude  (A)  en  tijd  (T).  Hele  lijnen 
representeren  T=20  min  dan  wel  F=0.1  Hz,  onderbroken  lijnen  T=60  min  of 
F=0.3  Hz  en  stippellijnen  T=120  min  of  F=0.5  Hz. 
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3  MODEL  VOLGENS  GRIFFIN 
3.1  Theorie 


Om  het  percentage  zeezieken  te  voorspellen  gebruikt  Griffin  de  term  “Vomiting  Incidence” 
(Vf).  Met  deze  VI  wordt  echter  precies  hetzelfde  bedoeld  als  met  de  term  MSI.  Om  een 
duidelijk  onderscheid  aan  te  brengen  tussen  de  resultaten  van  de  verschillende  modellen  zal 
hier  de  term  VI  expliciet  gekoppeld  worden  aan  het  model  van  Griffin.  Griffin  onderscheidt, 
evenals  McCauley  drie  grootheden  die  van  belang  zijn  voor  het  model:  de  mate  en  de 
frequentie(s)  van  de  beweging,  en  de  duur  daarvan.  Eventueel  kan  ook  de  gevoeligheid  van 
de  betrokkenen  als  parameter  in  zijn  model  worden  meegenomen. 

Mate  van  beweging:  Hiervoor  neemt  Griffin  de  versnelling  (a  in  m/s^).  Deze  maat  is 
evenredig  met  de  kracht  die  het  lichaam  ondervindt  door  de  beweging. 


Frequentie(s)  van  de  beweging:  Het  is  bekend  dat  mensen  bij  bewegingen  trager  dan  0. 1  Hz 
en  sneller  dan  0.5  Hz  nauwelijks  of  niet  bewegingsziek  worden.  De  grootste  gevoeligheid 
treedt  op  bij  bewegingen  met  een  frequentie  /  iets  onder  de  0.2  Hz.  Griffin  houdt  met  d^e 
frequentieafhankelijkheid  rekening  door  de  mate  van  beweging  te  wegen  met  een  factor  w 
die  als  volgt  afhangt  van  de  frequentie: 


w 


f 

0.125 


w  =  1.0 


1^  =  0 


voor 

0.1 

<  /  <  0.125 

voor 

0.125 

<  /  <  0.25 

voor 

0.25 

<  /  <  0.5 

elders 

(7) 


Niet  de  versnelling  a,  maar  zijn  gewogen  pendant  a„  =  w(f)  •  a  dient  nu  verder  in  de 
berekeningen  betrokken  te  worden. 

Duur  van  de  beweging:  De  kans  dat  iemand  moet  overgeven  onder  invloed  van  een 
beweging,  neemt  (uiteraard)  toe  met  de  duur  T  van  die  beweging.  Binnen  de  door  Griffin 
onderzochte  perioden  die  varieren  van  20  minuten  tot  6  uur,  is  gebleken  dat  de  MSI  niet 
recht  evenredig  toeneemt  met  de  tijd,  maar  na  enkele  uren  afvlakt. 


Griffin  definieert  nu  allereerst  een  dosiswaarde,  de  “Motion  Sickness  Dose  Value”,  waarbij 
het  cumulatieve  effect  tot  uiting  komt: 


(T 


MSDV  = 


/■ 


a},dt 


(8) 


Door  het  kwadrateren  van  de  versnelling  wordt  voorkomen  dat  bij  het  bepalen  van  het 
cumulatief  effect  opvolgende  negatieve  en  positieve  versnellingen  elkaar  teniet  doen.  Het 
nemen  van  de  wortel  uit  de  integraal  (het  cumulatieve  effect)  zorgt  ^n  voor  een  uiteindelijke 
evenredigheid  tussen  MSDV  en  versnelling  fen  voor  een  verzadiging  in  de  afhankelijkheid 
met  de  tijdsduur  zoals  uit  experimenten  is  gebleken. 
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Voor  een  enkelvoudige  sinusoidale  beweging  kan  vergelijking  2  vereenvoudigd  worden  tot: 


MSDV  =  w  • 


= 


w  •  a„ 


(9) 


Op  basis  van  deze  dosiswaarde  wordt  de  VI  berekend  met  vergelijking  10.  De  VI  is  ook  bier 
weer  een  cumulatieve  maat:  bet  totale  percentage  van  een  groep  dat  binnen  de  gegeven  duur 
moet  overgeven,  dan  wel  de  kans  dat  een  individu  dit  doet. 

VI  =  K  ■  MSDV  (10) 


Gevoeligheid:  Met  de  waarde  K  kan  de  gevoeligbeid  van  een  betrokken  groep  aangegeven 
worden.  Gemiddeld  beeft  voor  een  normale  burgerpopulatie  K  volgens  Griffin  een  waarde 
van  1/3.  Griffin  schrijft  vrouwen  een  5/3  maal  grotere  gevoeligheid  toe  dan  mannen^.  Wordt 
bier  expliciet  mee  rekening  gehouden,  dan  dient  K  de  waarde  0.4  respectievelijk  0.25  te 
hebben. 


3.2  Model  implementatie 

Voor  sinusoidale  stimulatie  kan  de  VI  met  vergelijkingen  9  en  10  eenvoudig  berekend 
worden.  Is  bet  bewegingsprofiel  samengesteld  uit  meerdere  frequenties,  zoals  bij  reele 
scheepsbewegingen  veelal  bet  geval  is,  dan  moet  voor  iedere  frequentie-component  de 
amplitude  gewogen  worden  met  (7)  waarna  pas  de  integraal  (8)  berekend  mag  worden.  Dit  is 
werk  voor  computers. 

Gegeven  een  equidistant  bemonsterde  tijdreeks  van  versnellingen  (in  m/s*),  dan  kan  met  een 
FFT-algoritme  snel  en  efficient  per  frequentie-component  de  versnelling  berekend  worden. 
Per  frequentie  kan  dan  de  weging  volgens  (7)  worden  uitgevoerd.  Na  terugtransformatie 
(inverse  FFT)  is  dan  een  bewegingsprofiel  gegeven  waarin  de  versnellingen  gewogen  zijn 
met  bun  frequentie-afhankelijke  weegfector.  De  integraal  over  deze  aangepaste  tijdreeks  kan 
benaderd  worden  door  de  cumulatieve  amplitude  (in  bet  kwadraat)  maal  de  tijd  te  berekenen. 
Door  toepassing  van  de  FFT  wordt  echter  wel  een  beperking  aan  de  data  opgelegd,  namelijk 
dat  bet  aantal  data-punten  een  macht  van  twee  dient  te  zijn.  Wanneer  dit  niet  zo  is,  moeten 
er  dus  data-punten  buiten  beschouwing  blijven,  of  punten  worden  toegevoegd.  Dit  laatste  kan 
bet  eenvoudigst  door  de  data  met  “nullen”  (een  interval  met  versnelling  nul)  aan  te  vullen. 
Bij  toepassing  van  deze  FFT  op  een  beperkt  computergeheugen  ontstaat  een  tweede 
probleem  wanneer  de  te  analyseren  data-set  te  groot  is.  In  dit  laatste  geval  kan  een  equi- 
distante  sub-set  genomen  worden  door  sub-sampelen  (bet  trekken  van  een  selecte  steek- 
proef),  of  door  bet  botweg  weglaten  van  de  laatste  sampels  (duur  inkorten).  Door  bet  sub- 
sampelen  kunnen  hogere  frequenties  ten  onrechte  worden  weggelaten  en  kunnen  aliasing 
effecten  optreden  (oorspronkelijk  niet  aanwezige  frequenties  werken  wel  door  in  de  analyse). 
Door  bet  inkorten  van  de  duur  wordt  op  triviale  wijze  de  VI  verkleind. 


^  Volgens  de  ISO  2631/3  norm  zijn  vrouwen  slechts  5  %  gevoeliger  voor  zeeziekte  dan  mannen. 
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Een  FORTRAN-programma  voor  MS-DOS  dat  de  MSDV  en  VI  berekent  zoals  hier  be- 
schreven,  is  inmiddels  gerealiseerd.  Een  voorbeeld  van  de  vragen  en  resultaten  zoals  dit 
programma  die  stelt  en  weergeeft  staat  in  Bijiage  B.  Het  input-file  dient  hierbij  een  enkel- 
koloms  ASCII-file  te  zijn  met  equidistant  bemonsterde  data-punten  die  de  versnelling  in  m/s^ 
geven. 


3.3  Resultaten 

Wanneer  de  beweging  sinusoTdaal  is,  kan  met  vergelijking  9  de  VI  berekend  worden.  Zo  is 
analoog  aan  de  berekeningen  aan  het  McCauley-model  ook  hier  voor  ieder  interval  van  0.1 
^  F  ^  0.5  Hz,  0.3  ^  A  ^  5  m/s^  en  20  ^  T  ^  120  minuten  de  V7  berekend.  Deze 
resultaten  staan  weergegeven  in  Fig.  4. 


Fig.  4  VI  versus  frequentie  (F),  amplitude  (A)  en  tijd  (T).  Hele  lijnen  represen- 
teren  T=20  min  dan  wel  F=0.1  Hz,  onderbroken  lijnen  T=60  min  of  F=0.3 
Hz  en  stippellijnen  T  =  120  min  of  F=0.5  Hz. 
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4  DISCUSSIE  EN  AANBEVELINGEN 
4.1  Discussie 

Het  grote  verschil  tussen  het  McCauley-  en  het  Griffin-model  is  dat  McCauley  e.a.  de  MSI 
voorspellen  enkel  voor  sinusoidale  bewegingen,  terwijl  Griffin  dat  ook  doet  voor  complexe 
bewegingen.  Beide  modellen  gelden  verder  enkel  voor  verticale  bewegingen.  Het  McCauley- 
model  is  geverifieerd  voor  frequenties  tussen  0.08  en  0.7  Hz,  het  Griffin-model  voor 
frequenties  tussen  0.1  en  0.5  Hz.  Voor  het  gebied  waarin  beide  modellen  geldig  zouden  zijn 
is  een  vergelijking  gemaakt,  waarvan  de  resultaten  in  Fig.  5  staan  weergegeven.  Hieruit  kan 
geconcludeerd  worden  dat  er  grote  verschillen  (tot  ongeveer  20%)  bestaan  tussen  beide 
modellen.  De  standaarddeviatie  van  het  absolute  verschil  tussen  beide  methodes  ligt  echter  in 
de  burnt  van  de  5%,  hetgeen  in  het  licht  van  de  grote  verschillen  tussen  mensen  onderling 
niet  dramatisch  is. 


Fig.  5  Verschillen  tussen  het  McCauley-  en  Griffin-model.  Hele  lijnen  repre- 
senteren  T=20  min  dan  wel  F=0.1  Hz,  onderbroken  lijnen  T=60  min  of 
F=0.3  Hz  en  stippellijnen  T=  120  min  of  F=0.5  Hz. 


Het  voordeel  van  het  Griffin-model  is,  zoals  genoemd,  dat  die  het  aantal  ze^ieken  voorspelt 
op  basis  van  bewegingen  die  meerdere  frequenties  bevatten.  Dit  is  vooral  voor  kleinere 
schepen  van  belang.  Het  Griffin-model  heeft  daarnaast  echter  een  aantal  niet  te  verwaarlozen 
nadelen.  HierbiJ  zij  op  voorhand  aangetekend  dat  dit  door  Griffin  zelf  al  wordt  onderkend. 
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Zijn  werk  is  vooral  gebaseerd  op  onderzoek  bij  passagiers  op  (grote)  veerboten.  Dez& 
passagiers  stellen  zich  i.h.a.  passief  op,  met  zicht  op  buiten,  tijdens  relatief  kort  durende 
overtochten.  Dlt  is  een  groot  verschil  met  scheepvarend  marinepersoneel  dat  zich  ook  actief 
gedraagt,  zelfe  benedendeks,  en  vaak  meerdere  dagen  aaneensluitend  op  zee  is,  maar 
waarschijnlijk  minder  gevoelig  is  voor  zeeziekte. 

Wanneer  voor  welk  model  gekozen  dient  te  worden  (Griffin  of  McCauley)  is  niet  precies  te 
zeggen.  De  beste  vuistregel  lijkt  mij  dat  wanneer  een  MS/-voorspelling  gewenst  is  voor  een 
enkelvoudig  sinusoidale  beweging,  of  een  duidelijk  unimodale  beweging  met  een  nauwe 
frequentie-band,  McCauley  de  voorkeur  verdient.  Dit  is  vooral  ingegeven  door  de  bevindin- 
gen  van  Colwell  (1994)  die  suggereert  dat  het  model  van  McCauley  e.a.  betere  (d.w.z.  meer 
met  de  door  Colwell  bemeten  werkelijkheid  overeenstemmende)  resultaten  geeft  dan  dat  van 
Griffin.  Op  basis  van  eigen  M5/-data  vergaard  op  zee  vindt  Colwell  dat  het  Griffin-model 
zijn  data  voorspelt  met  een  font  van  gemiddeld  4%  (sd=9%),  en  dat  het  McCauley-model 
dat  doet  met  gemiddeld  een  font  van  3%  (sd=7%).  Een  groter  nadeel  van  het  Griffin-model 
vindt  hij  echter  dat  voor  een  langere  duur  en  grote  amplitude  de  MSI  (=V7)  significant  wordt 
overschat  (zie  ook  Fig.  5).  Slechts  wanneer  het  bewegingsprofiel  complex  is  met  meer  dan 
ddn  duidelijk  overheersende  frequentie,  dient  het  Griffin-model  te  worden  toegepast. 

Wertheim  e.a.  (1995a)  hebben  al  aangegeven  dat  er  gegronde  redenen  zijn  om  te  twijfelen 
aan  een  model  waarin  alleen  de  domp-component  een  rol  speelt.  ZiJ  geven  een  voorbeeld 
van  33%  MSI  in  een  experiment  waar  Griffin  6%  voorspelt  (Wertheim  e.a.,  1995a).  Bij  een 
groep  als  niet-zeeziektegevoelig  bekend  staande  proefpersonen  vinden  zij  zelfe  50%  MSI 
waar  Griffin  3%  voorspelt  (Wertheim  e.a.,  1995b).  Zij  bevelen  dan  ook  aan  een  model  te 
maken  dat  naast  de  verticale  bewegingen  ook  het  slingeren  en  stampen,  en  indien  mogelijk 
ook  de  hoofdbewegingen  betrekt  bij  de  M5/-berekening. 

Een  laatst  te  noemen  tekort  van  beide  modellen  is  de  beperkte  geldigheidsduur  (tot  2  uur 
voor  McCauley,  20  minuten  tot  6  uur  voor  Griffin).  Bij  de  marine  zijn  juist  langer  durende 
vaarten  gebruikelijk,  en  het  is  bekend  dat  er  in  de  loop  van  dagen  gewenning  ontstaat  aan  de 
opgedrongen  beweging.  Wanneer  de  MSI  voor  de  opeenvolgende  dagen  opnieuw  wordt 
bepaald,  kan  deze  zelfe  tot  nul  reduceren  (Colwell,  1994). 


4.2  Aanbevelingen 

Ora  een  nadere  invulling  te  geven  aan  een  uitbreiding  van  een  model  ter  voorspelling  van  de 
MSI  met  slinger-,  stamp-  en  eventueel  hooffibewegingen,  is  het  volgende  vermeldenswaard. 
Uit  onderzoek  uitgevoerd  voor  de  Koninklijke  Luchtmacht  is  bij  de  programmagroep 
Evenwicht  &  Orientatie  bekend  dat  het  zgn.  Coriolis-efFect  verband  houdt  met  de  mate  van 
gevoeligheid  voor  luchtziekte  (Bles  e.a.,  1994).  Hierbij  maakt  een  persoon  terwijl  hij  draait 
(bijvoorbeeld  op  een  draaistoel)  een  kanteling  met  het  hoofd.  Dit  geeft  onaangename 
draaisensaties  van  bewegingen  die  in  werkelijkheid  niet  optreden,  hetgeen  het  Coriolis-effect 
genoemd  wordt.  Het  effect  is  soms  zo  hevig  dat  het  maken  van  de  hoofdbeweging  al  direct 
aanzet  tot  braken.  Bij  gecombineerde  scheeps-  en  hoofdbewegingen  zal  de  mate  van  dit 
Coriolis-effect  gering  zijn.  Het  is  echter  de  overtuiging  van  de  auteur  van  dit  rapport  dat 
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analoog  aan  vergelijking  8  ook  hierbij  een  cumulatie  optreedt,  die  uiteindelijk  tot  zeeziekte 
kan  leiden.  Het  ligt  dan  ook  voor  de  hand  MSI  mede  te  berekenen  op  basis  van  een  soort 
cumulatief  Coriolis-efFect. 

Zoals  ook  al  aangegeven  door  Weitheim  e.a.  (1995a)  kan  MSI  mogelijk  beter  voorspeld 
worden  m.b.v.  modelstudies  op  basis  van  kennis  van  het  evenwichtssysteem  (fundamentele 
aanpak),  dan  op  basis  van  kennis  van  MSI-scoros  bij  verschillende  bewegingen  alleen 
(empirisch/fenomenologisch).  Een  ruimtelijk  orientatiemodel  op  vestibulaire  grondslag  dat 
door  de  programmagroep  Evenwicht  &  Origntatie  voor  de  Koninklijke  Luchtmacht  is 
gemaakt  zou  daarbij  goed  als  basis  kunnen  dienen  (Bles  &  Bos,  1994;  Bos  &  Bles,  1995). 
Naast  een  “orientatie”-output  kan  dit  model  worden  uitgebreid  met  een  “onwelbevinden”- 
output.  Parameterkeuzes  worden  hierbij  dan  bepaald  door  bekende  MSI-scores  zodat  dit 
model  dat  van  McCauley  en  Griffin  overlapt. 
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BULAGE  A  Voorbeeld  MSI 

Hieronder  volgt  een  voorbeeld  van  de  output  van  het  programma  MSI,  versie  4-9-95,  zoals 
dat  ook  wordt  weergegeven  op  het  beeldscherm.  Antwoorden  (input)  volgen  na  een  het 
resultaat  (output)  volgt  na  de  “  = 


MSI  versie  4-9-95. 

Dit  programma  berekent  de  "Motion  Sickness  Incidence" 
volgens  McCauley  e.a.,  Human  Factors  Research,  Inc. 
Techn.  Rep.  1733-2  (1976). 

De  berekende  MSI  (%)  geldt  voor  verticale  sinusoidale 
bewegingen  met  een  RMS-ampl itude  tussen  0.0278  en  0.555 
G,  frequenties  tussen  0.083  en  0.7  Hz  en  een  duur  tot 
2  uur. 


RMS-Versne11ing  (G)  :  .2 
Frequentie  (Hz)  :  .3 
Duur  (minuten)  :  120 

Percentage  "spugers"  =  41.1 


EINDE  MSI 
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BULAGE  B  Voorbeeld  MSDV 

Hieronder  volgt  een  voorbeeld  van  de  output  van  bet  programma  MSDV,  versie  15-2-95, 
zoals  dat  ook  wordt  weergegeven  op  het  beeldscherm.  Antwoorden  (input)  volgen  na  een 
resultaten  (output)  volgen  na  een  “  =  ”,  de  rest  is  begeleidende  tekst. 


MSDV  versie  15-2-95. 

Dit  programma  berekent  de  "Motion  Sickness  Dose  Value" 
volgens  de  British  Standard  6841,  zoals  beschreven  in: 
Griffin  MJ  (1990)  Handbook  of  Human  Vibration,  Academic 
Press,  London. 

De  berekende  schatting  geldt  voor  verticale  bewegingen 
met  frequenties  tussen  0.1  en  0.5  Hz  en  een  duur  tussen 
de  20  minuten  en  6  uur. 

Input  file  met  z-versnell ingen  (m/s^)  :  test 

Sample- interval  in  sec.  :  1.0 

Inlezen  data  . . . 

Er  zijn  meer  data-punten  dan  toegestaan. 

Een  selecte  steekproef  nemen  (1), 

of  de  duur  inkorten  (2)  :  1 

Inlezen  data  . . . 

Er  zijn  9216  data-punten, 
waarvan  4608  gebruikt  met  een 
duur  van  2h  33min  36s. 

#  Gewogen  frequenties  =  4096. 

Frequentie  stapgrootte  =  .610E-4  Hz. 

Signaal  frequentie  <  .250E+0  Hz. 

[^0|^0|r]0p 

Motion  Sickness  Dose  Value  (MSDV)  =  2.6 

Percentage  "spugers"  =  .9 


EINDE  MSDV 


VERZENDLUST 


1.  Directeur  M&P  DO 


2.  Directie  Wetenschappelijk  Onderzoek  en  Ontwikkeling  Defensie 

Hoofd  Wetenschappelijk  Onderzoek  KL 

3.  { 

Plv.  Hoofd  Wetenschappelijk  Onderzoek  KL 

4.  Hoofd  Wetenschappelijk  Onderzoek  KLu 

Hoofd  Wetenschappelijk  Onderzoek  KM 

5.  { 

Plv.  Hoofd  Wetenschappelijk  Onderzoek  KM 

6,  7,  8.  Hoofd  van  het  Wetensch.  en  Techn.  Doc.-  en  Inform. 

Centrum  voor  de  Krijgsmacht 

9.  Ir.  E.M.  Krikke,  Hoofd  sie  Hydromechanica,  Dienst  Materieel  Koninklijke  Marine 
Afd.  Scheepsbouw,  Bureau  Wetenschap.  Zaken,  Den  Haag 


Extra  exemplaren  van  dit  rapport  kunnen  worden  aangevraagd  door  tussen- 
komst  van  de  HWOs  of  de  DWOO. 


